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ABSTRACT 

Context. The global properties of the interstellar medium with processes such as infall and outflow of gas and a large scale circulation 
of matter and its consequences for star formation and chemical enrichment are important for the understanding of galaxy evolution. 
Aims. In this paper we studied the kinematics and morphology of the diff'use ionized gas (DIG) in the disk and in the halo of the 
star forming spiral galaxy NGC 4666 to derive information about its kinematical properties. Especially, we searched for inf ailing and 
outflowing ionized gas. 

Methods. We determined surface brightness, radial velocity, and velocity dispersion of the warm ionized gas via high spectral res- 
olution (R ^ 9000) Fabry-Perot interferometry. This allows the determination of the global velocity field and the detection of local 
deviations from this verlocity field. We calculated models of the DIG distribution and its kinematics for comparison with the measured 
data. In this way we determined fundamental parameters such as the inclination and the scale height of NGC 4666, and established 
the need for an additional gas component to fit our observed data. 

Results. We found individual areas, especially along the minor axis, with gas components reaching into the halo which we interpret 
as an outflowing component of the diffuse ionized gas. As the main result of our study, we were able to determine that the vertical 
structure of the DIG distribution in NGC 4666 is best modeled with two components of ionized gas, a thick and a thin disk with 0.8 kpc 
and 0.2 kpc scale height, respectively. Therefore, the enhanced star formation in NGC 4666 drives an outflow and also maintains a 
thick ionized gas layer reminiscent of the Reynold's layer in the Milky Way. 
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1. Introduction 

^ It is well known that stellar winds and supernovae explosions 

"^ are the processes that contaminate the interstellar medium (ISM) 

CN with heavy elements ( [Shields 2002 ). Furthermore, supernova 

t^^ explosions can trigger additional star formation because of the 

1^ pressure they provide to the ISM. On galactic scales stellar 

t:j- winds and supernova explosions can influence the galactic 

^^ evolution in an additional way. They locally heat and accelerate 

CO the ISM, thus creating so-called superbubbles (SBs) (Tomisakaj 

^"^ |et al.|198T ). Because of lower pressures in the vertical direction, 

^ these SBs expand faster in the vertical direction creating 

• ^-H structures, much like smoking chimneys, through which the gas 

/\^ is able to flow into the halo of the galaxy ( [Norman & Ikeuchl] 

Depending on the energy of this outflow, one can distinguish 
two scenarios. Either the gas has enough energy to escape 
the host galaxy, creating a galactic wind that enriches the 
intergalactic medium (IGM) ( |Heckman| |2QQ2t [Veilleux et al.| 
2005 ), or the gas is unable to escape the gravitational potential 
of the host galaxy. In the second case the chemical enrichment 
of the galaxy is delayed and a large gaseous halo can be created 
aro und the galaxy. 

i Shapiro & Field (1976) described the processes of the gas 
dynamics into the halo and the fall-back as a galactic fountain. 
In this model the ejected gas condens es to clouds in the halo 
and falls back onto the disk. Recently, Marinacci et al. ( 2011| ) 



studied the galactic fountain model by using two-dimensional 
hydrodynamical simulations of a cloud ejected from the disk 
into the hot coronal gas. They described the influence of radia- 
tive cooling on the structure of this cloud and detected in their 
simulations a deceleration of the ejected cloud as a function 
of the height above the plane, i.e., a vertical velocity gradient. 
This is also called a lagging halo. Lagging halos are observed 
NGC 5775 dHeald et al. 



in several galaxies e.g 

NGC 891 ([Fraternali et all2005l JHeald et al.|2006a 



( [Fraternali et al.|2 
] |Qosterloo et al.|^ 



20Q6bl) 



or 



Kamphuis 



let al.|2007|TOosterloo et al.|2007| ) 

The ejected gas in the halo exists in diflTerent phases. This 
paper is focused on the extraplanar diflTuse ionized gas (eDIG). 
This extraplanar component of the diflTuse ionized gas (DIG) 
was detected in several galaxies, for example in our Milky Way 
(Re ynolds II 1990| ). The study of the eDIG in other galaxies has 
been made in several galaxies in Ha (,Dettm^|1990[ |Rossa & 



|Dettmar|2000l|2003a|bt|Rand et al.|1990t|Rand|[l996J ). Detailed 
kinematic studies of the halo gas were made in several dwarf 
galaxies ( ,Martin et al.||2002t |van Eymeren et al.||2009a|b| ) and 
continued in this paper for the SABc galaxy NGC 4666. 

Here we present a study of the kinematics of the diflTuse ion- 
ized gas in the disk as well as in the halo of NGC 4666 using 
Fabry-Perot interferometry. In this context we will determine the 
possible consequences of the strong starburst. This includes the 
detection of outflows, described by pahlem et al.| ( pW7| ), and the 
physical properties of an eDIG layer (thick disk). 
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Table 1. General parameters of NGC 4666. 



Parameters (Unit) 



Value 



optical center: 

a (J2000.0)^ 

6 (J2000.0)^ 

D [Mpc]^ 

Helioc. Radial velocity (km s~^)^ 

Galaxy type^ 

Size (Maj. x Min. Diameter) CY 

Inclination (°)^ 

Position angle {°y 

Max. rotation velocity (km s~^)^ 



12M5'"08.^6 

-00°27'43" 

18.225 

1529 

SABc 

4.6x1.3 

69.6 

42 

192.9 + 2.3 



Data from NASA/IPAC Extragalactic Database (NED) 

(http://ned.ipac.caltech.edu/). 

Data from the Lyon/Meudon Extragalactic Database (LEDA) 

( http ://leda.univ-ly onl .fr/ ) . 

Data from .Springob et al.| ( [2007| ). 



The gal axy NGC 4666 supp osedly has a high star 
: (SFR) ( [Walter et al.p 2004') of T.OMoyr"! (Dah 
^ ' - ' ^11 ■ ' ' ' ' 

er 



formation 



Dahlem et al. 



1991) . It has been well studied at several wavelengths, e.g., m 
Hot and R-band (Lehnert & Heckman 1995 , 1996a b ), radio con- 
tinuum and soft X-ra y (iDahlem et al.|1997||1998] ), and 21 cm and 
CO (Walter et al. 200^^ 

Garcia (1993) determined NGC 4666 to be a member of a 
small group of galaxies. Besides NGC 4666 this group is com- 
posed of NGC 4632 and NGC 4668. In this context it is currently 
believed that the strong st arburst is t riggered by ongoing far-field 
gravitational interactions ^ Walter e t al. 2004). Table [T] shows the 
general parameters of the galaxy. 

This paper is structured as follows: in Sect. 2 we describe 
the observations and the data reduction and in Sect. 3 follows 
the results of our analysis of the galaxy. In Sect. 4 our Yia obser- 
vations are compared to other observations, especially in other 
wavelengths, and in Sect. 5 we summarize our results. 

2. Observations and data reduction 

For our observations of NGC 4666 we used the high spectral res- 
olution (R^9000) Marseille Fabry-Perot interferometer (FPI) 
on the 1.93 m telescope at the Observatoire de Haute-Provence 
(OHP) in Franc e. The F PI was equipped with a photon count- 
ing system (Ga ch et al.| |2002) which provides a field of view 
(FOV) of 5.4 X 5.4 arcmin^ on 512 x 512 pixels^. This results in 
a pixel size of 0.63'' . Furthermore, we used a narrow band Ha fil- 
ter with a width of 10. 10 A and a maximum transmission of 70% 
at 6597.7 A. The systemic velocity of the galaxy is 1529 km s"^ 
(see Table[l]) which redshifts the Ha line to 6596.3 A, close to the 
maximum transmission of the chosen Ha filter. The free spectral 
range of the FPI is about 8.34 A or 379 km s"^ . This results, using 
a 32-channel mode, in a sampling step of 0.26 A or 1 1 .83 km s"\ 
respectively. By measuring the FWHM of a neon line, emitted by 
a neon calibration lamp, the finesse could be determined to 11.5 
and the spectral resolution to 0.73 A or 33.0 km s"\ which is less 
than 3 channels. 

Because of reflections within the FPI, a ghost image of the 
galaxy is produced. This ghost is mirrored in the center of the 
FOV. To prevent an overlap of the real signal and the signal of 
the ghost we positioned the galaxy about V off'set from the cen- 
ter of the FOV. This off'set is small enough that the galaxy does 



not leave the FOV, but is large enough to ensure a clear sep- 
aration between the ghost and the galaxy, thus preventing any 
contamination. 

The observations were obtained in two nights in early 
February 2011. In these two nights the galaxy was observed in 
50 and 47 cycles, respectively. A failure in the cooling system in- 
creased noise and created bad cycles, which were then removed, 
leaving us with 50 and 44 cycles for each night. To identify pos- 
sible extended halo structures we moved the telescope for the 
second observation about 1.7' west in order to position the ghost 
image once east and once west of the galaxy. With an integration 
time of 10 s per channel in one cycle we reached a total integra- 
tion of 30,080 s (^ 8.36 h) or 940 s (^ 15.7 min.) per channel. 

The data were reduced with the IDL program FPreduc 
( jDaigle et al.|2006[|Epinat|2009 |. Both cubes were reduced in- 
dependently. In order to correct for fluctuations in the telescope 
positioning, the position of every 10 s exposure was corrected by 
matching the positions of several stars in the FOV. Furthermore, 
the sky lines were subtracted and a heliocentric velocity correc- 
tion was applied. With the koords routine of the karma package 
(| Gooch| 1996) both cubes were provided with a coordinate sys- 
tem by matching four stars in the FOV to their counterparts in a 
Digitized Sky Survey (DSS) image. 

After the reduction with FPreduc the continuum was still 
visible in the data. Therefore, the continuum was subtracted us- 
ing the continuum image which was produced during the data 
reduction. Once both cubes were fully reduced they were moved 
to the same pixel coordinates and combined. 

The observed velocities of NGC 4666 cover a larger range 
than the available free spectral range. Therefore, some areas with 
extreme velocities appear in channels that are ofl'set from their 
true velocities by the spectral range. These velocities are mea- 
sured from the previous or next order of the FPI. In order to cor- 
rect this, the cube was extended with 16 channels, 10 channels 
in front and 6 channels behind the normal velocity range. Such 
parts of emission which contained velocities of the previous or- 
der were added in the channels behind the normal range and such 
parts of emission which contained velocities of the next order 
were added in front of the normal range. Furthermore, these re- 
gions were blanked within the normal range. 

The subsequ ent analysis of the data cub e was done with the 
GIPSY package ( |van der Hulst et al.||1992| ). Although GIPSY is 
mainly used for the analysis of interferometric observations of 
H I line emission, it has many functions that can also be applied 
to FPI observations of Ha emission. To optimize the S/N the 
cube was spatially smoothed by using the task SMOOTH, from 
3.6", which is equal to the average seeing, to a size of 5.0". 
This value of 5.0" was determined to be the best compromise 
between increasing the S/N in the outer parts of the galaxy and a 
good spatial resolution. For the following we define the number 
of complete pixels within a block of 5.0" x 5.0" as one resolution 
element, which corresponds to 7 x 7 pixels. 

3. Analysis 

In our analysis we focus on five diff'erent areas. Initially, we de- 
termine the sensitivity of our observation, then we verify the 
statement of (Wal ter et al.| |2004) that NGC 4666 is a starburst 
galaxy by determining the star formation rate. The third step is 
the velocity analysis of our data. From this analysis we obtain a 
rotation curve for NGC 4666 and a good impression of the kine- 
matics. In order to understand how much of the gas is in the disk 
and how much is in the halo, we use this kinematical informa- 
tion to construct a model for the disk. Finally, we compare our 
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disk model to the data to understand structures in the halo of the 
galaxy. 



3.1. Sensitivity of our observation and verification of the star 
formation rate 

To determine the sensitivity of our data we added all chan- 
nels to generate an integrated Ha map which we will re- 
fer to as the intensity map from here on. With a total 
Ua flux of 4.7 X 10"^^ erg s"^ cm"^ of NGC 4666, determined 



by Moustakas & Kennicutt| ( |20Q6| ) with drift-scanning spec- 
troscopy, we determined the 3cr sensitivity of our observations 
to be 5.6 X 10~^^ erg s"^ cm"^ arcsec"^. Another common unit to 
describe the intensity of emission is the so-called emission mea- 
sure (EM). It is given by 



EM 



=/"^ 



dr. 



(1) 



where n^ is the electron density and r the length along the line-of- 
sight. For a direct calculation we could use the formula ( [Dettmar 
|1992| ) 



EM = 5 X 10^^ F Ha/ ^cm-^pc. 



(2) 



with the measured flux Fna in ergcm"^sec"^ and D the col- 
lecting area in n''. This results in an emission measure of 
EM = 2.8 cm"^ pc for the aforementioned sensitivity. 

NGC 4666 is des cribed as a su perwind galaxy with a high 
SFR of 7.0MQyr~^ ([Dahlem et al.^1997 ) or even as a starburst 
galaxy ( Walter et al.||2QQ4 ). This value should also be verified 
for our discussion. For Ha we took the aforementioned flux of 
4.7 X 10"^^ erg s"^ cm"^ a nd for comparison the flux density in 
FIR at 60yum of 37.34 Jy ( |Soifer et al.| 19891 ). The Ua flux is in- 
fluenced by extinction. To obtain the extinction, we used t he HjS 



flux of 0.90 054 X 10 ^^ergs ^ cm ^ determined by Kennicutt 



etal. 
wEic 



( 2009 ), and calculated a B aimer decrement of ^ = 



5.22 



1 IS above the expected ratio of 2.86 for a temperature of 
10,0 00 K ([Osterb rock 1989). Taking into account the reddening 
law ( |Fitzpatrick| 1 999 ), the formula for the color excess E(B-V) 
is given by [Xiao et al.[ ( 2012j as 



E(B-V)= 1.99 X log 



HalHp 
2.86 ' 



(3) 



This results in a color excess of E(B-V) = 0.52 which 
corresponds to A(V) = 1.61, with A(V) = 3.1xE(B-V), 
and A(B) = 2.13. With A(Ha) = 0.538 xA(B) the extinc- 
tion of the Ha line can be calculated as A(HQf) = 1.15. 
Therefore, the total Ha flux is increased by a factor of 2.88 to 
13.51 X 10"^^ erg s"^ cm"^. For the additional calculations, we 
used the mean distance of 18.2 Mpc, provided in NEEQ whereas 



Dahlem et al. (1997) used 26.3 Mpc. The difl'erent values given 



by NED are mostly obtained by the TuUy-Fisher relation and 
are approximately 19 Mpc. The final value of NED we used is 
the average of all these single values. In contrast, the value of 
[Dahlem et al.| ( |1997| ) is based on the redshift with a correction 
for the Virgo infall. Since both methods are not very accurate 
we prefer to take the average from NED. Using the formulas of 
|Kennicutt| ( [T998l ), we obtained an SFR of 

SFR = 7.9 X lO'^^LHa ^"^^ = 4.24 Mo yr"^ (4) 

ergs"^ 



NASA Extragalactic Database,|http : //ned . ipac . caltech . edu/ 



for Hof and 

SFR = 4.5 X 10-^^4^ 



FIR 



Mpyr- 
ergs"^ 



= 2.2Moyr" 



(5) 



for far-infrared (FIR). 

Using the distance of 26.3 Mpc ( [Dahlem et"aLl|1997| ) the 
SFRs would increase by a factor of about 2. This results in a 
good agreeme nt between our Ha SF R and value of 7.0 Mq yr"^ 
determined by [Dahlem et al. ( [1997 ). The remaining diff'erences 
between the SFR in Ha and FIR can easily be explained by 
galaxy specific c haracteristi cs ([Kennicutt_1998j . In this context 
[Kennicutt[(l998[ ) referenced JBuat & Xu[ ( [i996[ ) who determined 
a coeflScient of 8^? x 10"^"^ instead of the above-mentioned 



4.5x10 



-44 



3.2. Kinematic analysis 

Before starting the further analysis of the data cube the above- 
mentioned ghost images of the galaxy below and above the disk 
were removed. For the further analysis the GIPSY task XGAUFIT 
was used. First a single Gaussian fit of every pixel with XGAUFIT 
was performed by using blocking filters for velocity and disper- 
sion/linewidtljj For the velocity only values within the veloc- 
ity range of the extended cube were accepted. For the linewidth 
the minimum value was set to the instrumental dispersion of 
14kms"\ which follows directly from the above-mentioned 
spectral resolution of the instrument. The maximal linewidth was 
set to 162 km s"\ which is equal to the dispersion of a Gaussian 
with a FWHM of the free spectral range. All these filters sup- 
press unphysical fitting results. The intensity map was masked 
at the 3(T level and subsequently used as a mask for the velocity 
field. As a last processing step a filter that blanked any isolated 
pixels was applied to the data. This results in a well-cleaned ve- 
locity field which can be seen in Fig.[T] All these steps were also 
done for the linewidth map which was produced with XGAUFIT, 
resulting in Fig. [2] (center). In both figures the black contours 
describe the intensity at 3cr, 6cr, I2cr, 24cr, etc. 

Because of our predominant interest in the vertical structure, 
the velocity map and the linewidth map were rotated so that the 
major axis of the galaxy is parallel to the x-axis. Moreover, the 
global coordinate system in right ascension and declination was 
replaced by a relative coordinate system that describes the radial 
and vertical off'set to the center of NGC 4666. 

These maps of NGC 4666 show a lopsidedness with a radial 
extent of the approaching part of about 110'' and a radial ex- 
tent of the receding part of about 125''. Furthermore, the veloc- 
ity map immediately shows several features beyond this lopsid- 
edness. When comparing the intensity contours in the velocity 
map, one can see the influence of the H ii regions to the mea- 
sured velocities, e.g., at the Hii regions at a radial/vertical ofl'set 
of-74"/16"and-30"/8". 

To determine the position angle and kinematical center and 
to obtain a rotation curve for NGC 4666, a tilted ring model was 
fitted by means of ;^^ minimization with the GIPSY task ROTCUR 
on the velocity map. Because of the lopsidedness of NGC 4666, 
only the inner 105" could be fitted well. The width of the rings 
was set to 5", our spatial resolution after smoothing. In a first run 
the position of the center of the galaxy, the rotation velocity, the 
systemic velocity, the inclination, and the position angle were 
all fitted in order to obtain initial estimates. Step-by- step single 
parameters were fixed using the average value over all the rings 



^ We define the linewidth as the dispersion of the Gaussian fit of the 
data which is not corrected for the instrumental dispersion. 
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Fig. 1. Velocity map of NGC 4666 with coordinates in offsets from the center of the galaxy in arcseconds. The contours show the 
intensity map at 3cr, 6cr, 12cr, 24cr, etc. 
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Fig. 3. The rotation velocity vs. the radial distance from the 
center. The blue triangles describe the approaching half of the 
galaxy, the red triangles the receding half, and the black dots 
both halves in a common calculation with the sigma of the veloc- 
ity given by ROTCUR as error bars. The semi-transparent triangles 
in blue and red describe the corrected rotation of the approaching 
and the receding half for the final model, respectively. 



to improve the fit for the other parameters. In the last iteration 
only the rotation velocity and the inclination were not fixed. The 
results are shown in Figs. |3] and |4] for the rotation curve and the 
inclination, respectively. All the other parameters and the aver- 
age of the inclination are given in Table [2j 

Looking at Fig. [3] one can see an almost flat rotation curve 
with a maximal rotation velocity of about 190-195 km s"^ which 
agrees well with the value found in the literature given in Table[T] 
The position of the galactic center and the position angle also 
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Fig. 4. The inclination vs. the radial distance from the center. 
Symbols as in Fig.|3] 



correspond to the values found in the literature (see Table [T]). 
The value of the position angle shows this offset of about 1 80° 
because of the definition in ROTCUR; in ROTCUR the position an- 
gle is defined as the counterclockwise angle between the north 
direction and the major axis of the receding part of the galaxy. 
The systemic velocity, the inclination and the position angle are 
the average values over all the rings up to 105'' with the standard 
deviation as the error. 

In the next step, the veloci ty offset along the mi nor axis of 
up to 100 km s"^ described by 



Dahlem et al. 



( 1997t win be in- 



vestigated. To verify this in our data we determined the velocity 
for every resolution element along the minor axis. Thus, we ob- 
tained for every vertical offset y along the minor axis a Av which 
is defined as Av = V3; - vq, with V3; as the velocity at a given verti- 
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Fig. 2. The linewidth map of NGC 4666, resulting from a single Gaussian fit, with coordinates given in distances in arcseconds to 
the center of the galaxy. The contours show the intensity map at Scr, 6(T, 12(T, 24(T, etc. Furthermore, some profiles of chosen pixels 
are displayed to show that high linewidths are double Gaussian profiles in the selected areas. When there are two Gaussians, the 
single Gaussian profiles are displayed in green and the sum of both in red. 



p. Voigtlander et al.: The kinematics of the diffuse ionized gas in NGC 4666 



Table 2. Results of ROTCUR for NGC 4666. The values are the 
averages of the inner ±105''. 



Parameters (Unit) 



Value 



kinematical center: 

a (J2000.0) 

6 (J2000.0) 

Systemic velocity (km s~^) 

Inclination (°) 

Position angle (°) 



12M5'"08.^6.8 
-00°27'43.9" 
1506 + 4 
78.3 + 0.5 
222.6 + 0.8 
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Fig. 5. The velocity offset along the minor axis of NGC 4666. 
The red tria ngles show our data; the blue diamonds show the 
data of Dahl em et al. ( |1997| ) . For this comparison the sign of the 
vertical offset in the data of Dahle m et al.| (|1997) was switched 
because of the different convention of the position angle. 



cal offset and vq as the velocity of the central resolution element 
(1510 km s"^). In Fig.|5]we plot A v vs. the vertical offse t of our 
data in comparison to the data of Dahlem et al. ( [1997 1. To do 
this, we switched the sign of the vertical offset in the data of 
[Dahlem et al.| ( [1997| ) because of the 1 80° difference in the as- 
sumed PA. 

In direct comparison, the general shape of the velocities be- 
tween -20'' and -h20'' vertical offset agrees well with the one 
presented injDahlem et al.| ( [l997| ). The velocity variations in the 
[Dahlem et al.| ( |1997| ) data and in our data in this inner area can 



Dahlem et al, 



\\99l\ ). 



this work, /Ires = 4 A = 182kms~^ m 

To study further this velocity offset in the outer regions and 
to detect other anomalies in the kinematics of the galaxy, the 
linewidth of the Gaussian fit is a useful tool. Therefore, we have 
chosen six areas to describe features in this linewidth map in 
more detail. The profiles of these regions can be seen in Fig. |2] 
around the linewidth map. First, the general linewidth is about 
30-35 km s"\ shown by way of example at position 1. The pro- 
file can be well fitted by one Gaussian (red line). Where we have 
compact H ii regions (black contour lines) the dispersion is about 
30 km s"^ or lower. In this linewidth map we see a region of 
large linewidth at a vertical offset of about 30'' along the mi- 
nor axis (position 2). The profile at this position shows a double 



be explained by the influence of Hii regions. The differences in 
the outer part are probably a combined effect of reduced S/N and 
a difference in velocity resolution (/l^ es = Q-73 A = 33 km s ~^ in [ 



Gaussian structure with a smaller Gaussian with a central posi- 
tion at 1421 km s"^ and a bigger one at 15 17 km s"^ (green lines), 
which is in good agreement with our central velocity, and also 
in better agreement with Dahlem et al. |1997) for this vertical 
offset (see Fig. [5]). By fitting merely one Gaussian to the line 
profile our central position is shifted to a lower value due to this 
first component at 1421 km s"^ Without that component the ve- 
locity keeps roughly constant for a positive vertical offset. This 
would be expected for a normal flat rotating galaxy. The first 
component with a velocity of about 100 km s"^ below the sys- 
temic velocity can be interpreted as an outflowing component of 
the gas. 

There are several other conspicuous areas in the linewidth 
map. In the center there is an x- shaped region with a larger 
linewidth which also shows a double Gaussian structure (e.g., at 
position 3). The linewidth map shows a linewidth for this struc- 
ture up to 60-80 km s"^ Furthermore, there is a significant rise 
of the linewidth, up to 50 km s"\ in the outer part of the galaxy 
above 20" and below -20" vertical offset. In the half of the 
galaxy with a negative radial offset the linewidth rises even up 
to 130kms"\ e.g., at position 4. Within these regions of large 
linewidths is a small hook-like structure where the linewidth 
drops down to 30-40 km s"\ marked as position 5. A detailed 
view of this position shows a similar shape as position 4, but 
here the amplitude of the first component, at 1335 km s"\ is sig- 
nificantly lower than the second component at 1427 km s"^ and 
close to the noise level. Therefore, the fit was dominated by the 
second component which has a small linewidth. Moreover, there 
is a region similar to position 2 where the linewidth rises within 
a compact area, marked as position 6 in Fig.|2] Both areas, 5 and 
6, will be discussed in more detail later in comparison to another 
independent Yia observation. 

The negative baselines of some profiles are a result of 
over-subtracted sky lines and because the continuum subtraction 
was optimized for the galactic center. Because of the over- 
subtraction, the above-mentioned calculated sensitivity is just 
an upper limit, but it does not affect the further profile analysis 
because it is merely a shift along the intensity axis and does 
not influence the shape of the profile. To understand the origin 
of all these double Gaussians we constructed a set of models 
to compare with the data. We will describe these models in the 
following section. 



3.3. Models 

To explain the aforementioned features in the velocity map and 
especially in the linewidth map and to determine the corrected 
rotation curve, an accurate model of NGC 4666 is essential. To 
model the data cube we used the FORTRAN code described in 



jKamphuis et al.j ( [2007 [ ). In this code the emission is calculated 
for every position in an exponential disk while taking the line of 
sight velocities into account. 

The model takes the following characteristic parameters as 
input. The intrinsic shape of the ionized gaseous disk is charac- 
teristed by an exponential scale length and scale height, as well 
as a truncation radius. The velocity structure is determined from 
an input rotation curve and the given linewidth. Finally, the disk 
is orientated according to a given inclination. 

The following parameters were determined directly from the 
data. The linewidth was measured from the diffuse gas in Fig.|2] 
as 35kms"^ The initial rotation curve and inclination were 
taken from the output of ROTCUR. The truncation radius could 
be determined easily from the total flux map as 130", which cor- 
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Table 3. Parameters for the best-fit model. 



Parameters (Unit) 


Value 


Exponential scale length (kpc) 


3.0 


Exponential scale height (kpc) 


0.2/0.8 


Inclination (°) 


78 


Truncation radius (kpc) 


11.5 


Distance (Mpc) 


18.2 


General linewidth (kms~0 


35 



responds to 1 1.5 kpc assuming a distance of 18.2 Mpc (see Table 

1). 

For the scale length and the scale height we started with some 
intitial estimates tracing the exponential intensity decrease along 
the major and minor axis, respectively. However, because of the 
large number of H ii regions in the disk of NGC 4666, these val- 
ues are very uncertain and had to be determined more accurately 
through a visual comparison of the data and the models. 

As a first step in fitting the model to the data, we corrected 
the rotation curve and determined an accurate value for the scale 
length by comparing the PV-diagrams along the major axis of the 
data cube and the model cube to each other iteratively. This re- 
sulted in an exponential scale length of 3.0 kpc and in corrected 
rotation velocities for the inner 15'' (see Fig. [5]). This correc- 
tion is not unexpected as we did not perform an envelope tracing 
( |Sofue & Rubin|200T] ) to determine the initial rotation curve. 

Subsequently, we determined accurate values for the incli- 
nation and the exponential scale height by calculating a grid of 
diff'erent models with diff'erent values for both parameters. This 
grid ranges from a scale height of 0.1 kpc to 1.0 kpc and from 
an inclination of 74° to 82°. The resulting model cubes were 
analyzed in the same way as the data cube which resulted in a 
sample of radial velocity and linewidth maps, as seen in Figs.|6] 
andJTJ respectively. 

We compared the models to the data through visual inspec- 
tion of velocity fields and linewidth maps. This confirmed an 
inclination of 78° on the one hand and revealed a more complex 
structure of the scale height on the other hand. The inner part 
of the galaxy matched well with the model with a scale height 
of 0.2 kpc, whereas the outer part matched well with the model 
with a scale height of 0.8 kpc. Therefore, we tried a combina- 
tion of these two models adding the thicker disk (z = 0.8 kpc) to 
the thin disk (z = 0.2 kpc). We found that a flux ratio between 
the thin disk and the thick disk of 2.5 matched the data best. A 
comparison of this final model with the data is shown in Fig. [S] 
for the velocity and linewidth map and confirms the quality of 
the calculated model. The parameters used in the final model are 
summarized in Table [3] 

We also investigated models with a dust layer and a lag- 
ging halo to study the possible influences of such components. 
However, these models did not show any indication of the pres- 
ence of a strong dust layer or a lagging halo. 

In Fig. 8] the white contours describe the model and the black 
contours the data. In the velocity map (Fig. [8] top) one can 
see again ofl'sets in velocity along the minor axis in the data. 
Otherwise, there is good agreement between model and data. All 
the residuals can be explained by the influence of H ii regions 
except the upper-right quadrant, especially at a vertical ofl'set of 
20''. The velocity in the data there is about 50 km s"^ higher than 
expected from the model. In comparison with our linewidth map 
(Fig. [2]), one can see a small rise of the linewidth; however, it 
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Fig. 9. The PV-diagram of NGC 4666 along the major axis with 
the data in black and the final model in white contours. 



remains in good agreement with our model. A possible explana- 
tion of this velocity off'set could be a warp. 

In the linewidth map of the final model (Fig. [8] bottom) one 
can see an x-shaped structure in the central region with a maxi- 
mum linewidth of 45 km s"^ whereas the linewidths in the outer 
parts show values of up to 60 km s"^ Compared to our data, the 
general structure is in good agreement. Nevertheless, there are 
two important diff'erences. First, in our data the linewidth in the 
center rises to 75kms"^ and is therefore 30kms"^ higher than 
in the model. Second, the rising linewidth in the outer part is 
well described by the model for a positive radial off'set; for a 
negative radial off'set the linewidth rises to 120kms"\ which is 
not described by the model. Increasing the inclination of the thin 
disk or thick disk does not solve this problem. We conclude that 
we are able to reproduce an increasing linewidth in an x-shaped 
structure in the center of the galaxy, but we are not able to reach 
the same strength. 

To study the quality of our model Fig. [9] shows the PV- 
diagram along the major axis. The chosen scale length of 
3.0 kpc agrees well with the data, but the assumed dispersion of 
35 km s"^ and the flat rotation velocity of 195 km s"^ also agrees 
well. 

As a result we conclude that the regions of large linewidth in 
our data, marked by way of example with 3 and 4 in Fig.[2j are 
partly explained as projection efl'ects of a thick and a thin disk 
in our data. For a complete explanation a detailed analysis and a 
comparison with other observations are necessary. 

3.4. Halo features 

To describe the halo features seen in Fig. [2] in more detail, a 
deeper view of the emission in the galaxy is necessary. For this 
we overlaid the Ha image from Dahlem et al.| ( |1997 ) with our 
velocity image in contours in Fig.[TO] To guide the eye, the black 
line representing a circle tilted at 78° with a diameter of 130" is 
drawn as a rough approximation of the disk size. 
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Fig. 6. Velocity maps of models with white contours in steps of 50 km s ^ around the systemic velocity of 1506 km s ^ at different 
inclination (y-axis) and scale heights (x-axis). The black contours describe the data with the same velocity steps. 



p. Voigtlander et al.: The kinematics of the diffuse ionized gas in NGC 4666 



E 



■D 

c 



- I I I I I I I I I I | >H I I I I I I I I - 

- ^'^ - 

- I i ^ 

- « ? » ^ 

L -tut : 



f 



MM IIII J. III MM 



MM MMJ.M I MM 



MM MMJ.M I MM 



■M 1 1 M 1 1 1 li^i 1 1 M 1 1 1 M- 



- 1 1 1 1 1 1 1 1 1 1 | j» 1 1 1 1 1 1 1 1 1 - 
^ t i ^ 






vt 



MM MMJ.M I MM 






MM MM J.M I MM 



■M 1 1 M 1 1 1 W\ 1 1 M 1 1 1 M- 



- 1 1 1 1 1 1 1 1 1 1 1 , 1 - 1 1 1 1 1 1 1 1 1 - 

- ^'^ - 

^ I i ^ 






lllllllll j.l. llllllll 

- { ^^ - 



wf 



lllllllll j.l. llllllll 

- i' '\ - 



-M 1 1 M 1 1 1 M 1 1 1 M 1 1 1 M- 



- 1 1 1 1 1 1 1 1 1 1 | j» i 1 1 1 1 1 1 1 1 - 
«* • 

f i - 



llll|llll j.U II|llll| 

- { ^.. - 

I i 



MM | MMjJ.M I| MM | 

ft S^^J^P^ 



MM | MMjJ.M I| MM | 

- { ^.. - 
! I 



M't 



vi 




■M 1 1 M 1 1 1 1'^l^i 1 1 M 1 1 1 M- 



- 1 I I I I I I I I I I y l I I I I I I I I I - 

- i' ^ - 

- 2 g - 

- • ^ J». — 

■mm. : 



llll|llll ].U II|llll 

- ^' ^.. - 

■mm. 



llll|llll ].U II|llll 

- ^'^ - 



~ { 



MM MMJ.M I MM 



-M 1 1 M 1 1 1 M 1 1 1 M 1 1 1 M- 



Q. 



Q. 

d ^-> 

II '-' 

N 53 



in 
it 
o 



■D 
03 

a: 



Q. 
i-l 

d 



O LD 



O ID 



O LD 



i = 74° 



i = 76° i = 78° i = 80° 

Vertical offset (arcsec) 



i = 82° 



Fig. 7. Linewidth maps of models with white contours at 45 km s ^ 55 km s \ 65 km s \ and 75 km s ^ at different inclination 
(y-axis) and scale heights (x-axis). The black contours describe the data with the same linewidth steps. 
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Fig. 8. The final velocity model (top) and the final linewidth model (bottom). The white contours show the model and the black 
contours the data. In the velocity map the contours are in steps of 50 km s"^ around the systemic velocity of 1506 km s"^ In the 
linewidth map the contours are at 45 km s"\ 55 km s"\ 65 km s"\ and 75 km s"^ 



First of all, we note that because of the lopsidedness of the 
galaxy the approximation of the disk with a simple ellipse fails 
on the left side. Furthermore, the above-mentioned influence of 
the H II regions in the velocity field is visible at several points, 
but more interesting are the structures below -25'' and above 
-f25'' vertical off'set. In the Ha image one can see extraplanar 
filaments reaching deep into the halo. In comparison to our ve- 
locity field one can see that our measurements go into these halo 
structures but are not as deep as the image by [Dahlem et al. 
( 1997| ). Nevertheless, a small protrusion in the velocity field at 
a radial and vertical off'set of 30'' also corresponds to an outgo- 
ing filament. The reasons that our observations do not show such 
structures in the halo are on the one hand the cut-off* due to the 
ghost images which would influence the structure in the outer 
halo and, on the other hand the necessity to have a significant 
signal above the noise to realize a Gaussian fitting. 

The structure of the galaxy in the lower-left quadrant is ex- 
tremely interesting. In all other quadrants the ellipse describes 
the shape of the galactic disk fairly well, but not there. Between 
the ellipse and the real disk is a distance of about 10" which 
corresponds to 0.9 kpc. A possible explanation is the lopsided- 
ness and, therefore, also a smaller size of the galactic height of 
the left part of the galaxy. However, in the linewidth map (Fig. 



[2]) this region shows a strong increase of the linewidth that can- 
not be explained with our models. Furthermore, in Fig. [TO] sev- 
eral filamentary structures and knots of higher emission like in 
the hook-like structure close to position 5 are visible. Therefore, 
these structures could represent outflowing extraplanar gas. 

The behavior of the velocity along the minor axis, especially 
the double Gaussian structure at region 2, can be understood bet- 
ter when taking into account the filamentary structures above 
-f25" vertical off'set that can be seen in the Ha image. We in- 
terpret this gas as part of an outflow. In region 6 no H ii region 
is visible, but strong emission of DIG up to a vertical ofl'set of 
-50". Therefore, an overlap of the normal rotation velocity and 
the velocity of the outflowing gas can be expected there. 



4. Comparison with observations in other 
wavelengths 

For the final interpretation a comparison with observations in 
other wavelengths is necessary. 

pahlem et al.| ( |1997|[T998] ) detected extraplanar soft X-ray 
emission and radio synchrotron emission from the halo and in- 
terpreted the se as further evidences for a galactic outflow. Walter] 
et al.| ( |2004l ) detected H i tidal streams connecting to NGC 4668. 
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Fig. 10. The Ha from |Dahlem et al.| ( 199 7) overlaid with the velocity contours of our data in steps of 25 km s ^ around the systemic 
velocity of 1506 km s~^ and a rough approximation of the disk size by tilting a circle with a diameter of 130'' to the inclination of 
the galaxy (78°). 



Furthermore, they assumed that the high star formation rate is 
triggered by gravitational interaction of NGC 4666 with the 
other galaxies in the group. For the Hi Walter et al.|(|2004|) de- 



termined an inclination for the inner 2' of about 70° whereas our 
determination of the inclination amounts to 78°. This results in 
a velocity offset of about 7 k ms~^ Taking this into account the 
maximal rotation velocity in Walter et al.| ([2004) is 173 km s"^ 
This is in good agreement to the value of the inner arcminute 
in our analysis of the ionized hydrogen and still in agreement 
to our maximal rotation velocity of 195 km s"^ when taking into 
account the errors. 



Although optical observations, including our own observa- 
tion, do not show any evidence of a strong interaction with other 
galaxies, such as a tidal streams, NGC 4666 cannot be seen as an 
isolated galaxy. Especially the lopsidedness of NGC 4666 could 
be a hint for such a gravitational interaction with other galax- 
ies. It is unlikely that the regular thick gaseous disk found in 
this study would be caused by environmental effects. But finally 
one has to take into account also the surrounding galaxies for a 
complete interpretation. However, the data for such an analysis 
is currently not available. 



5. Summary and conclusions 

We observed NGC 4666, an edge-on galaxy with a high star for- 
mation rate, with a Fabry-Perot interferometer to get line of sight 
velocity information of the ionized hydrogen in the disk and in 
the halo. Based on these data we calculated the rotation veloci- 
ties, which result in a flat rotation curve with a rotation velocity 
of about 195 km s"^ 

A comparison of the velocity along the minor axis, which 
should be constant, with the results of [Dahlem et al.| ( |1997| ) 
yields a good agreement within the inner ±20'' of vertical off- 
set. The differences at higher vertical offsets can be attributed 
to additional velocity comp onents, wh i ch are present in our data 
and not in the data of ,Dahlem et al.| ( |1997| ) as a result of our 
better spectral resolution. Especially for a positive vertical offset 
of more than 20'', a second velocity component with a velocity 
offset of about -100 km s"^ is confirmed, which we interpret as a 
strong outflow. 

To interpret the velocity and the linewidth of the fitted pro- 
files in an appropriate way we created models with different pa- 
rameters like inclination and scale height, based on the deter- 
mined rotation velocity. As result of comparing these different 
models with the data, a combination of two models, one with a 
thin disk (z = 0.2 kpc) and one with a thick disk (z = 0.8 kpc). 
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gave the best fit. A first important result is that the double 
Gaussian x-shaped structure in the center of the linewidth map 
as well as the double Gaussian structure in the outer regions at 
negative radial off'sets are results of projection eff'ects. The exact 
shape and the linewidth diff'ers on small scales between data and 
model, but because of the large number of H ii regions embed- 
ded in the diffxise ionized gas of the disk this is not surprising. 
We conclude as the main result of this paper the existence of a 
thin and a thick disk of diff'use ionized gas in NGC 4666, similar 
to the one in N GC 891 (|D ettmar"1990t |Rand et al.|[T9^ or in 
our own Galaxy ( [Reynold s 1990). Finally, we can conclude that 
modeling the linewidth map off'ers a proper possibility of detect- 
ing a thin and a thick disk of ionized gas in highly inclined, but 
not edge-on, galaxies. 

Furthermore, we found extraplanar structures, e.g., at a radial 
off'set of 25'' in both vertical directions, consistent with previous 
Ha observations. Clearly, the star formation in NGC 4666 main- 
tains a diff'use ionized thick disk and drives also a large scale 
outflow. If there is really a galactic superwind, cannot be con- 
cluded with our data. It is worth noting here that the outflow 
appears driven by a high star forming disk and not a nuclear 
starburst. This makes NGC 4666 more similar to z = 2 turbulent 
disks (e.g., Genzel et al.|2011 ) than merger driven nuclear star- 
bursts (e.g., |Heckman et al.|199Q| ). 
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